Persistent current in superconducting nanorings 
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The superconductivity in very thin rings is suppressed by quantum phase slips. As a result the 
amplitude of the persistent current oscillations with flux becomes exponentially small, and their 
shape changes from sawtooth to a sinusoidal one. We reduce the problem of low-energy properties 
of a superconducting nanoring to that of a quantum particle in a sinusoidal potential and show that 
the dependence of the current on the flux belongs to a one-parameter family of functions obtained 
by solving the respective Schrodinger equation with twisted boundary conditions. 

PACS numbers: 74.50.+r, 73.23.Ra, 73.63. Nm 



X 



The properties of superconducting grains change dra- 
matically when their size shrinks to a few nanometers. 
This phenomenon has recently generated a lot of in- 
terest Jj|. In particular, it has been established both 
theoretically and experimentally that the superconduc- 
tivity cannot be observed if the quantum level spacing 
in a nanoparticle exceeds the superconducting gap A. 
The recent success in manufacturing of superconduct- 
ing nanowires [0 has raised similar questions regard- 
ing the superconductivity in one-dimensional objects. It 
has been demonstrated experimentally || that as the 
nanowire becomes thiner, the superconducting transition 
in it disappears, and a finite resistance is observed at low 
temperatures. The suppression of the superconductivity 
in thin wires was attributed || to the destruction of the 
phase coherence by quantum phase slips 0, || . 

To understand the superconducting properties of thin 
wires one has to keep in mind that in an infinite one- 
dimensional conductor all the electronic states are local- 
ized. The localization length A in a wire of cross-section 
A can be estimated as A ~ kpAl, where k,F is the Fermi 
wavevector and I is the mean free path. If the length L 
of the wire is shorter than A, it can be viewed as a metal 
grain which becomes a good superconductor if the gap 
A exceeds the quantum level spacing Ep /kpAL. On the 
other hand, a wire of length L » A cannot be viewed 
as a good conductor, and its superconducting properties 
may be affected by localization. 

If the attractive interactions between the electrons are 
weak, so that the bulk superconducting gap A is small 
compared to the level spacing 5\ ~ Ep/kpAX of a piece 
of the wire of length A, the superconductivity is sup- 
pressed on a microscopic scale. In the more interesting 
regime of A 3> 6\, the superconducting gap is not af- 
fected by the localization pL The important question 
in this regime is that of phase coherence between the dif- 
ferent parts of the wire. Experimentally this issue can be 
studied by measuring the persistent current in a nanowire 
ring as a function of the magnetic flux piercing it. The 
magnitude of the persistent current oscillations as a func- 
tion of the flux is a direct measure of the superconducting 



phase coherence throughout the wire. 

In this paper we study the dependence of the persis- 
tent current in a superconducting nanoring as a function 
of the flux and the size of the ring. If the wire is rel- 
atively thick, the electrons move freely between differ- 
ent parts of the wire, and therefore the superconducting 
phase 4> is a well defined classical variable. At low tem- 
peratures T <C A there are no quasiparticle excitations 
in the nanoring. The allowed states of the system differ 
by the phase change ip = <fi(L) — 4>(0) accumulated over 
the circumference of the ring. At a given value $ of the 
magnetic flux through the ring <p can assume the values 
ip rn = 27r<E> / <!>q + 2irm, where m is an arbitrary integer, 
and $o = nfrc/e is the superconducting flux quantum. 
The energies of these states are given by 
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Here n s is the density of superconducting electrons, m* 
is the electron mass. The dependences of the energy lev- 
els ([!]) on the flux are shown by solid line in Fig. 0(a). 
The persistent current can be found as a derivative of the 
ground state energy I = cdE/dQ and shows the charac- 
teristic sawtooth behavior, Fig. |l|(b). 

The above picture fails in thin wires, where the fluc- 
tuations of the superconducting order parameter cannot 
be neglected. The most important type of the fluctua- 
tions is the quantum phase slip (J, [j| that changes the 
phase 4> at a point x by ±2ir. We will show that the 
effect of the rare phase slips on the system reduces to 
quantum transitions between the levels ([!]) . Such transi- 
tions are most important near the degeneracy points at 
half-integer values of $/§o and result in the small level 
splitting shown by dashed line in Fig. |](a). We will also 
show that the multiple quantum phase slips in longer and 
thiner wires eventually lead to a crossover from sawtooth 
to sinusoidal behavior of the persistent current, see dot- 
ted line in Fig. |l](b). 

The theory of quantum phase slips in nanowires [Q, |^] is 
rather complex. In particular, the most important quan- 
titative parameter describing the phase slip — its action — 
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FIG. 1: (a) The energy of a nanoring as a function of the flux 
$ through it. In the limit of relatively large cross-section A 
the energy shows the classical behavior ([l]) shown by solid line. 
Quantum phase slips result in level splitting shown in dashed 
line, and then lead to sinusoidal dependence of the ground 
state energy on the flux (dotted line), (b) The persistent 
current in the nanoring, /($) = cdE/dQ. In the classical 
limit the current shows sawtooth behavior shown by solid 
line. As the wire becomes thiner, the sawtooth is rounded 
and eventually transforms to a sinusoidal oscillation. 



can be found analytically only up to an unknown numeri- 
cal coefficient. On the other hand, the effect of the phase 
slips on the ground state of the nanoring can be studied 
using a much simpler model of a chain of coupled su- 
perconducting grains considered below. Similarly to the 
nanowire, the chain of grains will have quantum phase 
slips that affect the persistent current in the same way. 
We will discuss the relation between the two models in 
more detail below. 

We consider a model of the chain of superconducting 
grains defined by the following imaginary-time action: 
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The grains are assumed to be connected by tunnel junc- 
tions with Josephson coupling energy Ej] the capaci- 
tance of each junction C gives rise to the charging energy 
Ec = (2e) 2 /C; the variable 6 n is the phase difference 
across the n-th junction; (3 = l/T. To model a closed 
chain pierced by flux <I> = (tp/2Tr)&o, one should impose 
an additional condition 
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on the phases n . Note that the action (|^) docs 
not include the charging energy terms due to the self- 
capacitance of the grains. These terms, characterized by 
charging energy Ec per grain, are known j(| to give rise 
to a Kosterlitz-Thouless type quantum phase transition 



in an infinite chain of Josephson junctions at Ej/Ec ~ 1. 
The use of model (0) is motivated by the fact that the 
electric field is very well screened in metal wires which 
corresponds to the condition Ec ^S> Ej, Ec- 

The ground state properties of the chain can be derived 
from the partition function Z — J e~ s D6 n in the limit 
[3 — > oo. In this paper we consider the case of strong 
coupling between the grains Ej 3> Ec- At Ec = 
the phases O n become classical variables, and the energy 
states of the chain can be found by minimizing the sum 
of the Josephson energy terms in the action (|J) with the 
constraint (0). At large number of contacts N >• 1 one 
finds 



E m = (f + 2irmY 
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in complete analogy with Eq. (|l|) for a metal ring. 

At finite small Ec the fluctuations of 6 n (t) should 
be taken into consideration. First, one can include 
the Gaussian fluctuations around the classical solutions 
On(t) — const. This is accomplished by expanding the 
cosine in the action (|^) up to quadratic terms in 9 n . The 
resulting correction to the ground state energy SE rn = 
(N — l)-\/EjEc accounts for the zero-point oscillations of 
the particles in the minima of the cosine potential. This 
constant correction does not affect the shape of the oscil- 
lations of the persistent current / = (2e/K)dE/d<p. One 
can also include the anharmonic terms of the expansion 
of the cosine potential in the action (|J). The resulting 
corrections yield a small in Ec / Ej distortion of the saw- 
tooth dependence but do not significantly affect 
the discontinuities at half- integer <!>/<I>o. 

A more interesting fluctuation of the phases 6 n (t) in- 
volves an instanton (quantum phase slip), i.e., a trajec- 
tory that begins near one of the minima (^) of the poten- 
tial energy in action (|J) at t = and ends near another 
minimum at t = /3. For instance, a trajectory starting at 
0„(O) = tp/N and ending at 6 n (j3) = (<f- 2n)/N + 2-rt5 nk 
for arbitrary k in the interval 1 < k < N connects the 
minima (||) with m — and m = — 1. The shape of 
the instanton trajectory can be found by minimizing the 
classical action (^) with the above boundary conditions 
on 9 n {t). In the case of large number of junctions N ^S> 1 
the dominant contribution is due to the contact k where 
the phase slip occurs, and the contributions of the other 
contacts can be neglected. Then one obtains the usual 
result 
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where t' is the arbitrary moment in time where the phase 
slip is centered; we have assumed the limit of low tem- 
perature T <C \JEjEc- The action associated with this 
instanton trajectory is So = 8\/ Ejj Ec, 

The instantons account for the possibility of the sys- 
tem tunneling between the different minima @ of the 
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potential energy in action (g) . The effect of the instan- 
tons on the ground state energy can be accounted for by 
considering a tight-binding Hamiltonian defined as 



Hij) m = E m -)p m - Nv(ip m -i + VWi)- 



(6) 



Here ip m is the amplitude of the system at the state m 
with energy E rn given by Eq. (^). The hopping matrix 
element is exponentially small, 
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where the exponent coincides with the instanton action 
So, and the prefactor can be obtained by considering the 
problem in the case of a single junction (without the 
constraint (||)). 

It is important to note that for any k the set of phases 
6 n = (cp — 2tt)/N + 2n8 n k describes the same physical 
state of the chain. Thus the hopping matrix elements 
due to the instantons (|J) in all the N junctions must be 
summed up, resulting in the additional factor of N in the 
hopping term in the Hamiltonian (^). 

At Nv <C Ej/N the hopping term is small com- 
pared with the diagonal matrix elements in the Hamil- 
tonian @. Its effect is most significant when $/$o is 
half-integer, and the energy levels E m are degenerate, 
Fig. 0(a). In this regime hopping gives rise to the level 
repulsion shown by the dashed line in Fig. 0(a) and the 
respective rounding of the sawtooth in current, Fig. p]( 
The shape of the current steps is given by 
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where x — f - n = 27r($/<i>o - 1/2) < 1. 

In general hopping affects the spectrum of the Hamil- 
tonian (j(^) dramatically. To find the ground state en- 
ergy E(ip) it is more convenient to apply the Hamiltonian 
(||) to the wavefunction in "coordinate representation," 
^O) = E m ^ e ' (2m " lp/ ' ,)l - The resulting Schrodinger 
equation then takes the form of Mathieu equation: 

ip"(x) + (a- 2q cos 2x)ip(x) = 0, (9) 

where parameter a = 2NE/ir 2 Ej is proportional to the 
energy E, and q = N 2 v/tt 2 Ej. The phase <p enters the 
problem via a twisted boundary condition ip(x + n) = 
e~ lLf 'ip(x) . The regime of strong hopping in the Hamilto- 
nian (^|) corresponds to g > 1. In this case the depen- 
dence of the eigenvalue a on the phase tp is exponentially 
small, a = 



cos tp, Ref. @. For the 
persistent current / = (2e/h)dE/dtp we then find 



I = 16V2iV- r (£:, / u 3 ) 1/4 exp 
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It is important to note that the Schrodinger equation 
for the Hamiltonian (0) coincides with Eq. (||) at arbi- 
trary hopping strength q. Thus the problem of the flux 
dependence of the persistent current reduces to the solv- 
ing the well known Mathieu equation (|^) with the appro- 
priate twisted boundary conditions. This conclusion is 
our main result. The shape of the current oscillations is 
controlled by a single parameter q. As q changes from 
to oo, the shape crosses over from sawtooth to the sinu- 
soidal behavior. 

The Hamiltonian (^|) has been rigorously derived for 
the model (||) of a chain of coupled superconducting 
grains. Although nanorings are not described by this 
model, the essential low-energy physics of a supercon- 
ducting state being destroyed by quantum phase slips 
remains the same, and the Hamiltonian (^) and the re- 
spective Schrodinger equation (O) are still valid. 

To apply the Hamiltonian (0) to superconducting 
nanorings, one needs to relate the matrix elements E m 
and Nv to the parameters describing the nanowires, such 
as their length L, coherence and localization lengths £ 
and A, etc. The diagonal matrix elements E m are given 
by their classical values (jl]). The important difference 
between nanowires and chains of superconducting grains 
appears in evaluating the hopping matrix element Nv. In 
the absence of tunneling barriers the quantum phase slips 
must involve suppression of the superconducting gap A 
in the wire. The typical volume of space where the gap 
is suppressed is V ~ A£. Given the condensation energy 
per unit volume e ~ m*kp^ 2 and the typical correlation 
time r ~ 1/A, one can estimate the phase slip action 
in the wire So ~ eVt ~ y/A/6\ ~ A/£. A more rig- 
orous treatment |Q, || of the quantum phase-slips gives 
rise to the same result, but unfortunately the numeri- 
cal coefficient cannot be obtained analytically. Taking 
into account the obvious analogy between the number of 
Josephson contacts N and the length L of the wire, we 
conclude Nv cx Le~ So , with S — A/£. 

Another feature of realistic nanowires not included in 
our model of a chain of Josephson junctions is the meso- 
scopic fluctuations introduced by the disorder. For in- 
stance, since the instanton action depends on the mean 
free path, the fluctuations of the random potential can 
lead to the fluctuations of So. This effect should lead 
to mesoscopic fluctuations of the persistent current; the 
magnitude of the fluctuations is expected to be small in 
relatively thick wires. 

A more subtle effect of the random potential affects 
the relative phases of instantons in different parts of the 
wire. It can be explored within our model by adding 
gates to the grains and applying random gate voltages. 
It has been shown recently || [| that the gate voltages 
give rise to Aharonov-Casher interference effects for the 
instantons in different junctions, and under certain con- 
ditions can suppress the effects of the phase slips. 

To account for the effect of the random gate potential 
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on the persistent current, we present the electrostatic 
charging energy of the chain as 

H ° = E \[ c ^UiQn - q n )(Qi - a), (ii) 

nl 

where [C _1 ]„; is the matrix of inverse capacitances of 
the system, Q n is the charge at the n-th grain. Param- 
eters q n are determined by the gate voltages and have 
the meaning of the values of the charges Q n at which the 
electrostatic energy is minimized. 

Using the fact that the grain charge Q n and phase <j> n 
are conjugate variables and introducing 9 n = (j) n — (f> n +i + 
(fi/N, one rederives the imaginary-time action (||) with an 
additional term 



5S = -i^2,Pn 



Pn = qi + <?2 + • • • + qn 



(12) 



accounting for the gate voltages. Since the instanton (|^) 
corresponds to a trajectory where the phase in contact 
k changes by 2ir, the respective contribution to the hop- 
ping matrix element in the Hamiltonian m) has an ad- 



ditional phase e 



i27rp n 



Therefore for the total hopping 



matrix element we obtain vJ2 n e l2vPn instead of Nv, in 
agreement with Ref. |[ If the random gate voltages are 
sufficiently strong, so that q n ~ 1, the average of the 
hopping matrix element in the Hamiltonian (^|) vanishes: 
( v J2 n e l27rp ") = 0. Its typical values are of the order of 
the mean square fluctuation, thus decreasing from Nv in 
the regular case to ~ y/N v in the case of random charges. 
In the limit of large N, the same square-root dependence 
on the length of the chain will result even from weak dis- 
order. This suppression of the effect of phase slips results 
in a higher value of the persistent current, which can be 
estimated by replacing v — > v/y/N in Eqs. (0) and (10). 



The most significant change in the properties of the per- 
sistent current caused by the random gate potential is 
that instead of In / oc — N we find a slower dependence 
In I cx —N 3 / 4 on the length of the chain at large N. 

It is worth noting that the above result cannot be ob- 
tained from the action averaged over the random gate 
potentials q n . Indeed, one can easily see that (SS) = 0, 
thereby restoring our previous results (||) and (JlC|). A 
term similar to (|l2|) also appears in the theory of the per- 
sistent current in disordered wires. This term breaks the 
translational invariance and destroys the zero mode of the 
instanton action. Thus the conventional argument that 
the zero mode gives rise to the prefactor of the instan- 
ton contributions to the matrix elements proportional to 
the length of the system does not apply. Consequently 
we expect the total matrix element in the Hamiltonian 
(0) to behave as \fhe~ s °, leading to the persistent cur- 
rent In/ oc — Z 3 / 4 e~ s °/ 2 . Averaging of the action over 
the disorder would restore the zero mode and lead to the 
incorrect result In / cx —Le~ s °/ 2 . 

Our theoretical predictions can be tested in devices 
similar to those of Refs. ^| |3| Although we discussed the 



geometry of a nanoring of uniform thickness, all the con- 
clusions are applicable in the case of a straight nanowire 
with ends shorted by a bulk superconductor. In this ge- 
ometry only the quantum phase slips in the nanowire 
part should be considered, and thus the length of the 
system / is that of the nanowire only. The most inter- 
esting theoretical prediction to be tested in such devices 
is the dependence of the persistent current on the mag- 
netic flux The dependence /($) evolves from sawtooth 
to the sinusoidal one, Fig. g(b), as the nanowire becomes 
longer and/or thiner. We expect that the experimentally 
measured shape of the oscillations of /(<&) can be fitted 
to the function da(ip) / dtp, determined from the Mathicu 
equation (|^) . The only free parameter in such a fit is the 
effective phase slip strength q. 

Another type of systems to which our predictions 
should apply is the one-dimensional arrays of Joseph- 
son junctions. Recently it became possible to tune the 
Josephson coupling energy in such devices by external 
magnetic field jnj. As a result, one should be able to 
study the whole crossover from sawtooth to sinusoidal 
dependence of /(<&) in a single sample. 

In conclusion, we have studied the persistent current 
in nanoscale superconducting rings. As the wire becomes 
thiner, the quantum phase slips begin to suppress the 
current and change its flux dependence /($). The latter 
can be found by solving Mathieu equation (|^) with the 
appropriate boundary condition. The Aharonov-Casher 
interference suppresses the effects of the phase slips and 
gives rise to the unusual length dependence In / cx — Z 3 / 4 . 
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